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ABSTRACT 


A preliminary analysis of swell decay across the continental shelf was conducted 
using data from the SHOaling Waves EXperiment (SHOWEX). Six Datawell Directional 
Waverider buoys were deployed along a cross-shelf transect, offshore of Duck, North 
Carolina during an active hurricane season in the fall of 1999. Estimates of the frequency 
spectrum, E( f), mean propagation direction, Anegy( f), and directional spread, eT, 
were obtained from the auto-, co- and quadrature spectra of the horizontal and Genial 
orbital displacements measured by the buoys. Twelve cases were analyzed including 
eight cases dominated by remotely generated long period swells arriving from various 
directions. The last four cases were characterized by active wave generation by strong 
winds on the continental shelf, including the passages of iain Floyd and Irene. 
Comparisons of significant wave heights of swell observed at all buoys show a consistent 
strong decay across the shelf supporting the major role of bottom friction in swell 
transformation. A decline in wave heights near the coast was also observed during the 
sila of hurricanes, due to limited fetch and enhanced bottom friction in shallow water. 
Estimates of @meqn Show the expected refraction of waves to shore-normal directions on 


the inner shelf. Estimates of 6g generally do not vary much across the shelf. However, a 


sharp increase of Gg across the inner shelf was observed for directionally narrow swells 
from Hurricane Gert, possibly due to wave scattering from bottom irregularities. In all 


Cases, Og was at a minimum at the spectral peak frequency. 
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I. INTRODUCTION 


Swell characteristics transform dramatically between the open ocean and the 
beach. When the water depth becomes comparable to the swell wavelength, the bottom 
starts to influence swell propagation by causing changes 1n the wave phase and group 
speeds. | 

As swells begin to shoal the group speed initially increases, and it follows from 
the conservation of energy that the wave amplitude decreases by as much as 10% in 
intermediate water depths (kh ~1, sad k is the wave number and / is the water depth). 
In shallower water, the group speed decreases with decreasing depth, causing an increase 
in wave amplitudes near the coast. 

The phase speed changes associated with depth variations cause refraction 
analogous to Snell’s law in geometric optics. Monochromatic wave trains propagating - 
over a shelf with straight and parallel depth contours, refract such that the crests become 
more parallel to the beach. In situations with complex bathymetry, refraction can cause 
large wave height variations. Wave propagation directions generally refract towards 
shallower regions, typically causing a convergence region of larger waves near shoals, 
headlands and ridges, and divergence regions of smaller waves in bays and near 
submarine canyons (e.g., Munk and Traylor, 1947). Significant refraction effects occur 
as waves travel over a large submerged topographic feature, but also when waves move 
through a large shallow coastal region with smooth undulating eden (e.g., Graber 


et al., 1990; O’Reilly and Guza, 1993). 











Shoaling and refraction are linear processes that are well understood theoretically. 
By applying an energy balance along ray trajectories, the evolution of wave spectra can 
be evaluated for arbitrary two-dimensional continental shelf topography (e.g., Munk and 
Arthur, 1951; Longuet-Higgins, 1957; Le Méhauté and Wang, 1982). 

Although swell can cross entire ocean basins with minor energy losses (Snodgrass 
et al., 1966), significant dissipation of swell energy takes place on a shallow continental 
shelf. There are several mechanisms through which the bottom can remove energy from 
the waves and reduce wave heights. Shemdin et al. (1980) evaluated the relative 
importance of bottom friction (usually most dominant in sandy bottoms with 0.1-0.4 mm 
grain diameter), percolation through the bottom (coarser sediments with diameter = 0.5 
mm), and wave-induced bottom motion (fine sediment). The sediment type, bottom 
microtopography, in particular, the presence of ripples, and local currents play important 
roles in the damping rates, all of which are often unknown. 

On the sandy continental shelf considered in this study, bottom friction is believed 
to be the dominant swell decay mechanism. Bottom friction is dependent on both sand 
ripple characteristics and the prevailing sea-state conditions at the site. Tolman (1994) 
reviewed the important feed back between ripple formation and wave energy dissipation. 
Orbital sand ripples, initiated when the orbital wave motion near the bottom is strong 
enough to move sediment, greatly enhance the bottom roughness (Grant and Madsen, 
1982). Tolman’s (1994) model of wave-generated bedforms and associated damping of 
swell across a continental shelf suggest that moderately energetic swell shoaling on a 


continental shelf can generate sand ripples that cause strong damping of the swell. 





Few measurements of the attenuation of swell in shallow water have been 
reported. The JONSWAP experiment examined the swell decay observed over a 160- 
km-long cross-shore transect of various instruments deployed off the coast of Germany 
(Hasselmann et al., 1973). The observed damping rates did not agree with generally 
accepted formulations of bottom friction and suggested that other attenuation 
mechanisms such as back-scattering of swell from small scale bottom irregularities may 
be important (Long, 1973). Young and Gorman (1995) conducted a 2-week-long 
experiment using a cross-shore transect of various types of instruments spanning the 
continental shelf on the southern coast of Australia. These measurements were used in 
conjunction with the spectral wave model WAM (The WAMDI Group, 1988) to estimate 
the contribution of bottom friction to the observed overall decrease in swell energy across 


the shelf. The results suggested that the strong attenuation of energetic southern ocean 


swells observed at shallow instrument sites is primarily caused by bottom friction with 


estimated drag coefficients that are much larger than the values used in wave prediction 
models. Forristall and Reece (1985) measured wave attenuation and bottom processes 
near the Mississippi Delta in the Sea Wave Attenuation Measurement Program 
(SWAMP). They suspected that the soft bottom near river deltas causes strong swell 
attenuation. The observed attenuation was a strong function of the deep-water significant 
wave height and a weak function of the wave frequency. Refraction and shoaling 
accounted for the changes in the directional spectra when the wave height was small. 
With larger wave heights (> 3 m), significant attenuation was observed, in particular for 


frequencies < 0.10 Hz. 


Herbers et al. (2000) studied the ecidiinaien of swell across the continental 
shelf off the coast of North Carolina. The results suggested that on the outer shelf energy 
was lost only by the lower frequencies of the spectrum, possibly because high-frequency 
components with relatively shorter wavelengths were unattenuated at the bottom. At the 
shallower sites, observed spectral levels were lower than predicted spectral levels at all 
frequencies, suggesting significant damping of all components of the spectrum. The 
apparent damping is strongest near the peak frequency sinising a slight broadening of the 
observed spectrum at the closest nearshore site. The inferred friction factor agreed 
roughly with the results from Young and Gorman’s study, and was much larger than 
those commonly used in wave prediction models (e.g., The WAMDI Group, 1988). 

New observations of the evolution of swell across a continental shelf are 
presented here. Six Datawell Directional Waverider buoys were deployed along a 100 
km cross-shelf transect off the coast of North Carolina for 3 months in 1999 as part of the 
SHOWEX (SHoaling Waves EXperiment). In contrast to previous experiments that 
provided estimates of swell height decay across the shelf, the detailed directional wave 
measurements collected in SHOWEX also allow for an analysis of the transformation of 
directional wave properties across the shelf. An unusually active hurricane season caused 
by La Nifia conditions provided numerous case studies of energetic long period swell. 
The field experiment is described in Chapter II, followed by a description of the data 
analysis procedures in Chapter III. Twelve case studies are presented in Chapter IV, 


followed by an analysis of observed cross-shelf variations in Chapter V. Chapter VI 


discusses the frequency-directional spectra, and the results are summarized in Chapter 


Vil. 
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Il. FIELD EXPERIMENT 


The field data used in this study were collected as part of the ONR SHOaling 
Waves EXperiment (SHOWEX) conducted offshore of the U.S. Army Corps of 
Engineers Field Research Facility (FRF) near Duck, North Carolina, from mid-September 
to mid-December, 1999. The site is located about midway between Cape Hatteras and 
Chesapeake Bay, along a coast consisting of a series of relatively straight barrier islands 
with sandy beaches that are fully exposed to the Atlantic Ocean. 

The dikes shelf is oriented approximately north-south, and widens from 
Cape Hatteras toward Chesapeake Bay. The continental shelf is about 100 km wide and 
only 20 — 50 m deep across most of the study site (Figure 1). Bathymetry is characterized 
by a gently sloping inner shelf (about 1:250), extending into a wide mid-shelf ~— with 
an even lower slope of about 1:2500 and large (O(5 km)) ridge-like features that have 
vertical scales of O(5 m). The continental shelf widens from Cape Hatteras, North 
Carolina, and curves to the northeast toward Nova Scotia and the Grand Banks. 

Six Datawell Directional Waverider buoys (X1 to X6) were deployed along a 
cross-shelf transect extending almost due east from the Duck beach to the shelf break 
(Figure 2). The buoys were deployed in water depths of 21m (X1, located 5.35 km from 
shore); 24 m (X2, 12.78 km); 26 m (X3, 22.17 km from shore); 33 m (X4, 37.74 km from 


shore); 39 m (X5, 58.22 km from shore); and 193 m (X6, 86.57 km from shore) (Table 


1). 











Buoy Depth Latitude Longitude Dates 
P| prom Shore a aa 
fom Lee 36° 13.627N | 75° 42.28’ 9/13/99 — 
ee ee Aan399 

24m | 36 13.62’N | 75° 36.77’ 9/11/99 — 
ae ee ee 

x3 22.2 km 26m | 36 12.277 N | 75° 29.917 Wi 9/10/99 — 
Fe ome bine pres 

a lekeeiel 09.50°N | 75° 19.53’ W} 9/10/99 — 
ee EES ian 


er 2 km 36 07.79" N | 75° 06.05’ W | 9/13/99 — 
10/24/99. 
10/28/99 — 
12/13/99 
86.62 km 197m | 36 04. 
9/21/99 
sein 


36° 04.977 N | 75° 47.50’ W | 10/15/99- 
12/5/99 

Table 1. Locations and time periods of Datawell Directional Waverider buoys 

(X1 to X6) deployments 


























Buoy X5 broke loose on October 24, and was redeployed on October 28. Buoy X6 broke 





loose on September 21, was redeployed on October 15, and lost again from December 5 
through the end of the experiment. All other buoys were operational throughout the 
experiment. 

| The Datawell Directional Waverider buoys are surface following directional wave 
buoys (about 1 m diameter) which measure horizontal and vertical orbital displacements 
at the sea surface from which wave height spectra and directional wave properties are 
determined. They contain a heave-pitch-roll sensor (Hippy-40), three axis fluxgate 
compass, two fixed ‘x’ and ‘y’ accelerometers, a temperature sensor and a 
microprocessor. The direction measurement is based on the translation principle, using 


the measured horizontal accelerations rather than surface slopes used in earlier pitch-and- 








roll type buoys. Therefore, the measurement is independent of buoy roll motions and 
thus, a small spherical buoy can be used. From the accelerations and tilts measured in the 
moving “buoy reference frame,’ the accelerations along the fixed horizontal, north and 
_ west axes are calculated. All three accelerations (vertical, north and west) are digitally 
integrated to get filtered displacements with a high frequency cut-off at 0.6 Hz. To save 
transmitting power, the real time data are compressed to vertical, north and west 
displacements. Every half hour, Fast Fourier Transforms of 8 series of 256 data points 
(200 sec) are used to form spectra and cross-spectra with 16 degrees of freedom. 
Frequency resolutions of 0.005 Hz in the range 0.025 to 0.1 Hz, and 0.01 Hz in the range 
0.1 to 0.58 Hz are used. For each frequency bandwidth, the on-board microprocessor 
computes energy density, mean direction, directional spread, and the normalized second 
harmonic of the directional distribution (discussed in Chapter III). The Directional 
Waverider buoy transmits the following fafonnntion: vertical, north, and west 
displacements, spectral density, processed directional acuiiee sea surface temperature 
(SST) and other system status with parity bits for error checking purposes. The buoy 
data, framed in 64 bit vectors, 1s transmitted via HF in the 27-40 MHz band continuously. 
The processed spectral and directional data is transmitted 8 times to allow for reliable 
data return even in marginal transmission conditions. 

A single point vertical mooring (Figure 3) ensures sufficient symmetrical 
horizontal buoy response even for small motions at low frequencies. The Directional 
Waverider is fitted with 5 kg chain ballast to provide stability when only a vertical 


mooring force is present (free floating or moored in shallow water). The low stiffness of 


the 30 m rubber cord allows the buoy to follow waves up to 40 m in height. High shock 
accelerations are absorbed by the fluid surrounding the accelerometer and the shockproof 
mounting of the accelerometer housing, electronics and batteries. 

The four buoys located closest to shore (X1 to X4) transmitted their respective 
data via an HF radio link to receivers mounted on a 50 m high tower at the FRF, while 
the outer two buoys (X5 and X6) that were out of HF range, used internal data loggers. 

Significant weather associated with La Nifia conditions played an important role 
in the formation of long period energetic swell that propagated across the continental 
shelf during the fall of 1999. Hurricane season peaked in September and October with 
Hurricanes Floyd, Gert, Irene and José (Figure 4). On September 17, the eye of 
_ Hurricane Floyd passed within 100 km of the buoy transect generating large seas on the 
continental shelf and causing massive flooding and destruction in eastern North Carolina. 
Hurricanes Gert and José remained offshore but radiated long period swell across the 
transect of buoys. Hurricane Irene made landfall across western Cuba and southeastern 
Florida and continued northeastward, following the Gulf Stream and remained offshore 
of the Carolinas. Further discussion on these storms is found in the case studies in 


Chapter IV. 
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HI. DATA ANALYSIS 


The directional properties of ocean surface waves are often described by the mean 
wave. direction, 9,,.., (7), and the directional spread, o,(/), as functions of frequency. 
These variables are closely related to the directional aianibation of wave energy, 

S(O; f): Onean(F ) is the mean of this distribution, defining a dominant wave propagation 


direction at each frequency, whereas o,(/) represents the standard deviation of the 
distribution and characterizes the half-width of the directional distribution. Estimates of 
6 en (f)and o,(/) are readily obtained from heave-pitch-and-roll or translational 
directional wave buoys (e.g., Longuet-Higgins et al., 1963; Mitsuyasu et al., 1975; Long, 
1980, Kuik et al., 1988). Longuet-Higgins et al. (1963) showed that the first four 
directional Fourier coefficients of S(@; f) can be derived from the auto-, co- and 
sadeaare spectra of the three basic heave and slope (or horizontal displacement) signals 
of the buoy. For a Datawell Directional Waverider buoy that measures the horizontal 
(x,y) and vertical (z) displacements of water particles at the surface, the Fourier moments 


can be expressed as: 


O0..(f) 


ee: sr ] 
{CA MICa( A+ Cy (A 


a(f)= [dO cos (6) S(6;,f) = 


1] 


0,,(f) 


b, y= [40 sin (6) S(O: = ——___~3-2______ 2 

(A= Ja0 sin) SOD = oe H+ @) 
_" pA) 

a,(f) =| 48 cos (28) S(O 1) == Cf) (3) 
ae 2C,,(f) 

b,(f)= [6 sin (26) S(O; f) =——~ 4 

(f) J sin (20) S(O; f) A+ (4) 


where C;j and Q;; denote the co- and quadrature spectra of components i and j 
respectively. _ 


For narrow S(@), @,,,,ando, can be approximated by the first-order moments, 


a, and B, (e.g., Kuik et al., 1988): 


[ao sin (0-6,,.., ) S(9) = 8, COS (@rrean) —% SiN (Aoean) = 9 (5) 


a2 = [dO (21-605 (8-Oyeu)1) S() = 21 (a, C05 (Open) +4 Si (Orex))- 6) 


For small 0-6... 








sin (0-0...) © (8-Oye = = (0-B yn yy 


2[1—cos (8-Bgan))* (8-Bge) = (8-8 yen) 


and thus, @,.,, and o, represent a mean value and standard deviation of the directional 


distribution of wave energy with an O(0—@,.,)° bias. 
A disadvantage of this formulation is the sensitivity of o, to directionally 


Opposing wave components. A small amount of wave energy traveling in a direction 


opposing the dominant swell system (e.g., reflection of waves from shore or a weak 
secondary swell peak) causes a large increase in theo, estimate that is not related to the 
beam width of the dominant swell (Herbers et al, 1999). 


Alternatively, 0,.,, and o, can be defined in terms of the second-order moments 


a, and b,(e.g., Herbers et al., 1999): 


[ao sin [2(0-6.,.. )] S(8) = stb, cos (26,.,.)—a, sin (26,,.,,)] =0 (7) 


0} = [dO 1 {{1-208 [2(0-Oyan)}} S18) = 311 -(@, £05 (20x) +8, Si0 (28g) (8) 


For small 0-8... 
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5 Sin [2(B-Bgan)]* (8-Oan == (8-Ore 


5 {1~C05{2(8-Bgan)D} © (8-Oyan) = (8-Open 


and therefore the second-order approximations of 6_,, and o, have larger biases than 
the first-order approximations, Eqs 5 and 6. However, Eqs 7 and 8 are much less 
sensitive than Eqs 5 and 6, to a weak directionally opposing wave system, because the 
double-angle argument, 2(0-0,,.,,. ), is the same for waves propagating in opposite 
directions. Unlike the first-order estimate of 0... (Eq 5), the second-order estimate has a 
180° ambiguity. 

In the present analysis, estimates of 0,.,, and o, were obtained from both the 
first- and second-order moments. Because the two methods have different errors, 
intercomparison of @_.,, and o, estimates based on a,, b, and a,, b, (discussed below) 
provides a useful check of the accuracy of the estimates. A magnetic correction of 
—10.75° W was applied to 6,.,, (7), and the resulting direction was neaicaae, 
nautical convention, that is the transformed @__,(/) is the direction the waves come 
from, measured counterclockwise from true north. In addition to @,,,, and o, estimates 


as functions of frequency, bulk estimates were computed for the frequency range 0.025 to 
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0.25 Hz. The estimates are based on Eqs 5-8 using energy-weighted average values of 


the moments, a,, 5, and a,, b,: 


f aE 4 


q, = 208M (9) 


E(f) df 


0.025 Hz 


0.25 Hz 


[ a NE) of 
b ; — 0.025 Hz (10) 


0.25 Hz 


> (11) 


0.25 Hz 


| b,, tal ee (12) 


where E(f) [= C_(/)] is the wave frequency spectrum. 
A preliminary analysis of the bulk wave and wind parameters during the entire 


- experiment was conducted to select specific events for further analysis. For each case 


15 


covering time periods of 12 or 3 hours, the spectra of the half-hour data records were 


ensemble averaged, and directional moments, 6@,..,(/) and o,(/), were computed based 
on the ensemble averages of the moments (a,, b,, a,, and b,) weighted by E(f). The 


resulting E(f), 0...,(/), and o,(f) have approximately 384 degrees of freedom. 
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IV. CASE STUDIES 


The times of the selected twelve case studies are listed in Table 2, together with 


the significant wave height: 


H, 4 | Bas, (13) 


the peak period, T,, (the wave period corresponding to the peak of the spectrum, E(/)), 
and the bulk @,,., and o, parameters (calculated using Eqs 5-12) all from the offshore 
buoy, and the wave height decay across the shelf quantified by the ratio H,(X1)/H,(X6). 

Eight of the cases were collected during light wind conditions on the shelf when 
the wave field was dominated by remotely generated long period (T,=12-17 s) swell. 
Since it takes several hours for the swell to traverse the shelf, 12-hour periods were used 
in the analysis with relatively steady significant wave heights, peak periods and mean 
propagation directions. 

Six of the swell cases are “easterly swell” events with mean offshore directions 
within +/- 45° of 090° T. The easterly swell cases include forerunners from Hurricane 
Floyd and arrivals from Hurricane Gert. Hurricane Gert was a category IV hurricane that 
was tracking to the northwest and slowly turned north-northeastward on September 21. 


Gert’s closest point of approach was about 1100 km to the southeast of the buoy transect 


17 








(Figure 4). Gert produced the longest period (T,= 15-17 s) and directionally narrowest 


(a, =20-24°) swells of the experiment. 


EASTERLY SWELL H, at X6 | T, at Omen at 5 at H,(X1) 
X6 H,(X6) 
09/14 00 EST to 09/14 12 EST | 2.90m | 13.335 — ~ T | 29. 0.57 
(FLOYD) 


09/19 12 EST to 09/20 00EST | 2.78m | 15.38s | 1243°T | 204° 0.47 
(GERT) 


09/20 21 EST to 09/21 09EST | 2.54m | 16.67s | 127.0°T | 244° 0.55 
GRD 

09/23 16 EST to 09/24 04 EST XS: XS: XS: X5: | (XIX): 
_ STS [atm [12508 | mer | ase [0s 


11/18 03 EST to 11/18 15 EST 1429s | 545°T | 365° | 0.66 
11/22 18 EST to 11/23 06 EST 97.0°T | 348° | 0.67 


LARGE OBLIQUE ANGLE foe. oF 
10/23 02 EST to 10/23 14 EST 170m | 12.50s | 163.4 T 
Moe) ee 


12/01 21 EST to 12/02 09 EST 0348°T | 312° | 050 
Sowwaes WAVES | 




















a OD EST to 09/16 12EST | 7. — m | 12.50s > OT 0.62 
| (FLOYD) 
10/18 00 EST to 10/18 03 EST | 430m | 10.00s | 161.3°T 0.62 






(IRENE) 
11/02 20 EST to 11/02 23 EST 10.00 s | 167.9°T 
| T1123 EST to 11/12 02 EST [358m | 9.09s | 0203°T 


Table 2. Twelve case studies classified as: “Easterly Swell’, “Large Oblique Angle 
Swell” and “Storm Waves” listed at the offshore buoy, X6, are: Significant wave height, 
H,; peak period, T,; bulk mean propagation direction, mean; and bulk directional spread, 
O,. Also listed is the wave height decay between X1 and X6. 






In the remaining two cases, swell arrived at very large angles to shore-normal | 
~ resulting in strong refraction effects. The two “large oblique angle swell” cases include 
southerly swell from Hurricane José and a northerly swell case. Hurricane José remained 
well offshore and followed a somewhat similar track to Hurricane Gert, traveling to the 


north-northwest until passing the eastern tip of Puerto Rico before recurving to the north- 
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northeast (Figure 4). José’s closest point of approach was 900 km of the buoy transect to 
the southeast. The southeasterly swell from José propagated sues the steep slope of the 
continental shelf located to the southeast of the buoy transect (Figure 1) near Cape 
Hatteras Point. In contrast, waves arriving from northerly directions travel a long 
distance over the continental shelf, which widens and curves to the northeast. 

In addition to the 8 swell cases, 4 storm cases with strong winds and active wave 
generation on the continental shelf, were selected for analysis. Because of the more 
rapidly varying conditions, 3-hour periods at the peak of the storm events were analyzed. 
These cases include Hurricanes Floyd iat Irene, the most significant storm events during 
Pn entire 3-month experiment. Their tracks passed within 100-200 km from the buoy 
transect and produced large waves (maximum H, of 4.3 — 7.6 m) on the continental shelf. 
Figure 5 displays the associated winds and H, of Hurricanes Floyd and Irene. The peak 
winds of 35 m/s for Floyd occurred around 2000 EST on September 16 (well after 
Hurricane Floyd made land fall) while the peak H,; at buoy X6 were observed about 2 
hours later (Figure 5a). Since Floyd passed to the west of the buoy transect (Figure 4), 
the strongest winds were observed closest to shore as indicated by the difference in 
maximum wind speeds between FRF and NOAA Buoy 44014 observations. The NOAA 
Buoy 44014 was located about 30 km northeast of the buoy transect (Figure 1). 
Hurricane Irene made landfall over southeastern Florida, but remained offshore of the 
Carolinas, and passed within 250 km of the buoy transect (Figure 4). As Irene passed to 


the southeast of the transect, the dominant waves continued to arrive from the south- 
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southeast even though the strong local winds (on the west side of the hurricane) were 


from the northwest (Figure 5b). 


The other two storm events ied were local storms (a cold frontal system on 
November 2 and a nor’easter on November 11-12) that produced comparably large waves 
on the shelf. In the Floyd, Irene and the November 2 storm events, when the dominant 
seas observed at X6 arrived from the south, strong refraction and fetch-limited conditions 
on the inner shelf were expected. The nor’easter case shows very different conditions of 
fetch-unlimited seas generated by northeasterly winds. 

The November 2 storm moved slowly across the study site with wind speed 
observations from the FRF at a maximum of 18 m/s at 1200 EST while the NOAA 
buoy’s peak wind speed occurred about 10 hours later, during the 3-hour period selected 
(Figure 6a). In contrast, the nor’easter storm quickly passed through the study area with 
the peak wind observations at the NOAA buoy occurring 3-4 hours later than at the FRF. 
Hs was at a maximum at buoy X6 about 2 hours after the peak wind observations from 
the FRF station (Figure 6b). 

For each case the transfer function, G( f/f), was examined between the horizontal 


and vertical displacements, to verify the buoy response: 


Ic C 
G(f)= Sartre (14) 


For linear waves G(f) = 1 / tanh( k(f) h ), where h is the water depth and the 


wavenumber, k, is given by the dispersion relation (27 f)* = gk tanh(kh). Hence, the 
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normalized transfer function G'(/) [= tanh (k(/) h) G(/)] is equal to 1 for linear 
, waves. Nonlinearities in both the wave field and buoy response will result in deviations 


of G’ from 1. Bulk values of G and G’ integrated over the spectra are expressed as: 








G= (15) 
[ And 
0.025 Hz 
| f tanb(k(f) 2) [C,.9) + C, Alef 
Gl = |0.025 He ——— ee ee (16) 
| Ad 


Figure 7a shows the observed ratio, G, for each of the twelve cases. The observed ratios 
are close to 1 at the deep water (buoy X6), but, for all cases, G generally increases from 
the outer shelf to the inner shelf. This trend is consistent with the transformation of the 
wave vial motion from deep (circular orbits, G = 1) to shallow (ellipsoidal orbits, G > 
1) water. The normalized ratios vary between 1 to 1.2 (Figure 7b), ile a slight but 
systematic positive bias that increases with distance from shore, possibly due to nonlinear 


effects. 
The accuracy of the directional moments, 6... ando,, , is verified in Figure 8 by 


2? ™ mean 


comparing the estimates based on the second-order approximations (a,, b,) for each of 
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the twelve cases listed in Table 2, with their respective estimates based on first-order 
approximations (a,, 5, ). The bulk @ ean Values obtained with the two different 
approximations are in excellent agreement (Figure 8a). However, the o, values are well 
correlated but estimates based on (a,, b,) values are generally about 15-25% smaller 
than estimates based on (a,, 5, ) (Figure 8b). These differences are qualitatively 
consistent with the larger negative bias expected for o, estimates based on the second- 
order moments discussed earlier. The estimates of 6,..,ando, used in Chapters V and 
VI for the preliminary analysis of the twelve cases are based on the less-biased first-order 


approximations (a,, 5, ). 
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V. CROSS-SHELF VARIATIONS 


Variations in wave properties across the continental! shelf can result from different 
attenuation mechanisms discussed in Chapter 1. oe friction, refraction of waves 
arriving at oblique angles, and fetch limitations of wave growth on the continental shelf 
all contribute to a decrease in wave height across the shelf, causing sheltered conditions 


near the coast. In this Chapter, a preliminary analysis of the observed variations in H,, 


8 


mean ? 


and o,, across the shelf is presented. 


A. EASTERLY SWELLS 

The cross-shelf variations for the six easterly swell cases (Table 2) are displayed 
in Figures 9 and 10. In all six cases, a large decrease in wave heights across the shelf is 
observed, with H,; at buoy X1 on the inner shelf that are between 33 and 53% smaller 
than H, values measured with buoy X6 on the outer shelf. These decreases in wave 
height are much larger than the variations resulting from shoaling and refraction (Herbers 
et al., 2000), suggesting that the strong attenuation of swell is caused by bottom friction. . 
The strongest ‘iscay is observed for the most energetic swells generated by Hurricanes 
Floyd (September 14) and Gert (September 19 and 20). Although the overall decay is 
comparable for these two events, the cross-shelf evolution is different. During Hurricane 
Floyd, strong: initial swell decay (~37% decrease in H,) is observed across the outer shelf 
between buoys X6 and X5, followed by weaker variations across the mid- and inner 
shelf. During Hurricane Gert, a dramatic decrease in H, (between 25 and 37%) is 


observed on the inner shelf between buoys X3 and X2, with weaker variations offshore of 
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buoy X3 and inshore of buoy X2. The remaining three cases show a more gradual wave 


decay across the shelf. 
The easterly swell cases are characterized by bulk 6... ranging from 050° to 135° 
T (Figure 10a). The coastline and continental shelf break are oriented at ~ 070° T and 


100° T, respectively, and thus, the observed @,,.,, correspond to small (< 45°) angles of 


incidence. Variations of 6,,,,, across the shelf are small, except for the more southerly 
Floyd and Gert cases that show a ~ 30° shift toward normal incidence. In most cases, the 
bulk o,, remained either fairly steady (between 35° and 28°) or show a slight narrowing 
across the shelf that is likely due to weak refraction effects. The cases associated with 
Hurricane Gert with narrow offshore spectra (o,, 20°) show a sharp increase of o, of 
30% across the inner shelf (buoy X3 to X1) (Figure 10b), possibly due to wave scattering 
from bottom irregularities. 
B. | LARGE OBLIQUE ANGLE SWELLS 

In the two oblique angle swell cases shown in Figure 11, swells approach the 
shelf break at near grazing (~ 65° relative to normal incidence) northeasterly (December 


11) and southeasterly (October 23) angles. In both cases the significant wave height, H; 


decreased by about 50% across the shelf (i.e., energy reduction of a factor of 4) (Figure 


lla). This decay is likely caused by a combination of refraction and bottom friction. 


In both oblique angle cases, refraction effects are evident in the gradual turning of 
0 ea, across the continental shelf toward smaller oblique angles (Figure 11b). However, 


the observed refraction of the southeasterly swell from Hurricane José is much stronger 
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(from 165° T at buoy X6 to 105° T at X1) than in the northeasterly swell case where the 
change is less significant with swells arriving (from 035° T at buoy X6 to 055°T at X1). 
These direction shifts are qualitatively consistent with refraction computations presented 
in Herbers et al. (2000, cf. their Figure 3d). 

The o, (Figure 11c) observed in the northeasterly case is fairly uniform across the 
shelf (about 30° to 40°). In contrast, the southeasterly swells are broader at the offshore 


buoys, X5 and X6 (o, ~ 50-55°) narrowing to 38° at the inner shelf buoy X1. The 


unusually large o, offshore does not appear to be related to the shelf topography (buoy 
X6 is in deep water) and suggest possible scattering effects of the Gulf Stream. 
C. STORM WAVES 

The cross-shelf evolution of seas generated by Hurricanes Floyd, Irene and the 
local November 2 storm are presented together in Figure 12. The dominant waves 
arrived from south-southeast directions in all three cases, indicating fetch-limited 
conditions on the shelf. The nor’easter storm is an unlimited fetch case with — 
arriving from unobstructed north-northeast directions. Therefore, the nor’ easter case is 

presented separately in Figure 13. 

In the hie eetiodeas cases, the observed wave heights decrease across the 
continental shelf by about 40-50% between buoys X6 and X1 (Figure 12a). These H, 
reductions on the inner shelf are comparable to the decay observed in swell conditions 
(cf. Figures 10, 12a), but the mechanisms may be quite different as waves are actively 
generated on the shelf. In all three cases, the 0,,,, eStimates show waves arriving from 


| south-southeasterly directions at the outer shelf buoys, X6 and XS, shifting by about 60- 
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65° to near shore-normal directions at buoy X1 (115-120° T). The large differences 
observed between the dominant wave and wind directions on the inner shelf (e.g., 


southerly winds and easterly waves during Floyd) can be explained by the extremely 


limited fetch for southerly winds and the strong refraction of waves propagating nearly 


parallel to the local depth contours (Figure 12b). In all three cases a, ~25-30° at the 
outer shelf buoys (X4-X6), but the estimates diverge on the inner shelf (Figure 12c). In 
Hurricane Irene o,, increases to 45° at buoy X1, whereas the o, values of Hurricane 
Floyd and the local November storm case decrease slightly to 20-25° at buoy X1. These 
differences are discussed in Chapter VI. Variations of Hs, @,..,, and o, observed in the 
nor’easter case are small (Figure 13). H, decreased by only 25% across the continental 
shelf. The observed @,,,,, are about 020°-025° T at buoys X6 and X5, shifting to 040°- 
045° T at buoys X4 to X1. The o, are nearly uniform across the buoy transect varying 


from 25°-30°. These weak variations indicate that refraction is of lesser importance in a 


nor’easter, and that the wave field is fully developed in equilibrium with the wind. 
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VI. FREQUENCY DIRECTIONAL SPECTRA 


The evolution of the feqieneydiecionl wave spectrum, S(@; f), across the 
continental shelf is examined in detail for the four storm wave events: Hurricanes Floyd, 
Irene, the November 2 storm, and the nor’ easter. 

Estimates of E(f), Oren (f), and o,(f) at the peak of Hurricane Floyd are 
shown in Figure 14. The frequency spectra, E(f), (Figure 14a) show a large decrease in 
spectral levels between offshore buoy X6 and inner shelf buoy X1 in the frequency range 
0.08 to 0.12 Hz that corresponds to the peak of the offshore spectrum. The spectral peak 
level decreases monotonically from the outer shelf to the inner shelf (about a factor of 5 | 
between buoys X6 and X1), even though the maximum winds of 35 m/s were observed at 
the FRF station close to the inner shelf buoys (Figure 5a). At both higher frequencies and 
lower frequencies, the spectrum is nearly uniform across the shelf. The decrease in wind 
sea peak energy on the inner shelf may be caused by bottom damping and/or limited 
fetch. 

The @...,(/) estimates show large cross-shelf variations caused by refraction 
(Figure 14b). At the outer shelf buoys, X5 and X6, southerly angles (180° T through 
200° T) are observed at all frequencies. Between buoys X5 and X1, a shift in direction is 
observed that is largest for the lower frequencies. This is qualitatively consistent with 
refraction that is eon for low frequency waves that feel the bottom over the entire 

shelf. For example, at buoy X1, the mean direction at high frequencies is about 190° T 


(close to the local wind direction), and decreases with decreasing frequency to 120° T 
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(near normal incidence) at the spectral peak. Topographic effects play an important role 
in the evolution of lower frequency waves, but have little effect on the high frequency 


waves that are in equilibrium (1.e., saturated) with the local wind. 


For all six buoys, o,(/) is narrowest at the spectral peak (15°-25°), and broadens 


at both higher and lower frequencies (Figure 14c). No systematic variations of o,(/) 
across the shelf are observed, but at most frequencies, the spreads are generally narrowest 
at the inner shelf buoys (X1-X3) and widest at the outer shelf buoys (X4-X6). 

Although Hurricane Irene passed to the southeast of the buoy transect, producing 
different conditions on the shelf compared to Hurricane Floyd, there are many similarities 
in the spectra. The spectral levels were not as high as during Floyd, but the variations of 
spectral levels across the shelf are similar with about a factor of 5 reduction between X6 
and X1 at the spectral peak and weak variations at lower and higher frequencies. 

At low abies the @,..,(/) estimates of Irene are similar to those observed 
during Floyd, showing angles of refraction from about 180° to 120° at inner shelf buoy 
X1 (Figure 15b). At high frequencies (> 0.15 Hz) a wide range of directions is observed. 
At X1, high-frequency waves arrive from the northeast while at X6, south-southeasterly 
angles are observed (Figure 15b). These differences can be explained based on the local 
wind observations. At the time of the analysis, Irene was just passing to the south of the 
buoy transect and the winds were rapidly shifting from south to north at the FRF station, 
while wind measured at the NOAA buoy remained in the southerly direction until much 


later in the day on October 18 (Figure 5b). The high frequency part of the spectrum 
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apparently responds rapidly to the sudden shift in wind direction on the inner shelf, 


resulting in a shift to northeasterly directions at buoys X1 and X2. 
The corresponding o,,(/) estimates show a large increase in spreading at high 


frequencies between the outer shelf (30-45°) and the inner shelf (55-70°) that is likely 

. caused by bimodal directional spectra on the inner shelf with new northeasterly seas and 
old southeasterly swells. At lower frequencies, the o,(/) estimates are similar to those 
observed in Floyd with smaller values at the spectral peak frequency and weak variations 
across the shelf. 

In the November 2 storm case, the energy spectrum decreases uniformly across 
the shelf (Figure 16a). At all frequencies, spectral levels at buoy X1 are almost an order 
of magnitude smaller than those observed at X6. This large deeeease is consistent with 
large variations of winds across the shelf. The outer shelf experienced much stronger 
southerly winds (19 m/s at the NOAA Buoy) than the inner shelf (9 m/s at the FRF) 

_ (Figure 6a). Estimates of 6,,..,(/) (Figure 16b) show a large shift from southwesterly _ 
directions on the outer shelf to southeasterly directions on the inner shelf. At buoy X6, 
high frequency components propagate in roughly the local southerly wind direction. But 
at buoy X1, the local wind (200° T, Figure 6a) and high frequency wave directions 
(135°T, Figure 16b) differ by 75°. These observations indicate that high frequency 
waves observed on the inner shelf are generated by stronger winds on the outer shelf and 
refracted over the sloping shelf bathymetry. At high frequencies, the o,(/) estimates 


show a systematic decrease across the shelf (Figure 16 c) with much narrower spread at 


buoy X1 (25°-30°) , compared with buoy X6 (45°-50°), that is qualitatively consistent 


pay 


| | 


with the expected directional narrowing of waves refracting from grazing angles toward 
normal incidence. 

The nor’easter case on November 11 shows weaker variation across the shelf 
compared to the previous cases. At high frequencies, spectral levels are uniform across 
the shelf, characteristic of a saturated spectral tail in equilibrium with the wind. At low 
frequencies (< 0.15 Hz) spectral levels decrease monotonically across the shelf, 
suggesting attenuation of long waves by bottom friction (Figure 17a). At the peak 


frequency (0.11 Hz), the spectral decay between buoys X6 and X1 is about a factor of 2. 
At all buoys, 6... (/) is within the range 15°-45°, close to the wind direction, except at 


very low frequencies (cf. Figures 17b and 6b). The o,(/) estimate is very similar at all 
six buoys, with smallest values of 20-25° at the spectral peak, broadening at lower and 
higher frequencies (Figure 17c). | 

The directional spread observed for all twelve cases is summarized in Figure 18. 
The common trend for all buoys is a minimum spread at the spectral peak (// f,= 1), and 
increasing spread for both lower (f/f, < 1) and higher (f/f, > 1) frequencies. These 
observations are consistent with deep water observations by Mitsuyasu et al. (1975) and 
Hasselmann et al. (1980). The range of spread values at f= f, (listed in Table 3) does not 
vary much between the outer and inner shelf buoys. At higher frequencies the observed 
directional spreads are more variable on the inner shelf (cf. buoys X1 and X6), reflecting 
the complex wind generation (i.e., limited fetch) and topographic effects (e.g., strong 
refraction along coast hugging waves, scattering effects) near the coast that can both 


increase or decrease the directional spread. 
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Table 3. Directional spread range for buoys X1 to X6 at /=1 in Figure 18. 
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IV. CONCLUSIONS 


Duritie the SHOWEX (SHOaling Waves EXperiment), a cross-shelf transect of 
six Datawell Directional Waverider buoys a deployed on the North Carolina 
continental shelf to investigate the transformation of wind-sea and swell spectra in 
shallow water. A wide range of wave conditions was observed during the three month 
deployment including several major hurricanes. A preliminary analysis of the buoy 
observations is presented here. Estimates of the significant wave height, H,, frequency | 
re E(f), mean direction, Omean( /), and directional spread, o9( f) were extracted 
from the.buoy measurements using standard techniques. The evolution of wave 
properties across the shelf is examined to evaluate the relative importance of bottom 
friction, refraction effects, and limited wind fetch. 

The analysis is focused on twelve case studies that include remotely generated 
swell from Hurricanes Floyd, Gert and J osé, energetic wind sea conditions during the 
passage of Floyd and Irene, and a typical nor’easter for this region. In all swell cases, a 
large (33-53 %) decrease in wave height across the shelf was observed that confirms the 
donna role of bottom friction in the transformation of swell across a wide continental 
shelf. These results _ consistent with earlier apaeeaidons at the same site (Herbers et 
al., 2000). 

Estimates of Omean( Ff) show the gradual refraction of waves toward shore-normal 


directions. For swells arriving at the shelf break at a large oblique southerly angle, Omean 
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shifts of 7 to 90° are observed on the inner shelf, enhancing the wave height decay 
across the shelf. 

During the passage of Hurricane Floyd, wave heights decreased from 7.8 m at the 
shelf break * 4.9 m on the inner shelf even though the strongest winds were observed 
near the coast. Possible mechanisms that contribute to the wave height decrease are 
enhanced bottom friction in shallow water, limited fetch on the inner shelf for southerly 
winds and strong esfacion Similar evolution was observed during Hurricane Irene. In 
contrast, weak wave heights decay (23% decay) between the offshore and inner shelf 
buoys) was observed during a nor’easter storm. Waves are more uniform in this case 
because the fetch is unlimited for northeasterly winds and refraction effects are weaker 
for waves arriving from northeasterly directions. 

For all twelve cases at all six buoys, the oe( f) estimates consistently show a 
minimum (25° to 30°) near the spectral peak frequency. Cross-shelf variations of o are 
generally sitll with the exception of the directionally narrowest, Hurricane Gert swell 
observations that show a dramatic broadening across the inner shelf possibly owing to 


bottom scattering effects. 
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APPENDIX. FIGURES 
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Figure 1. Locations of Datawell Directional Waverider buoys X1 to X6 across the 
continental shelf. | 
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Figure 2. Bottom profile along the cross-shelf transect. The locations of the Datawell 
Directional Waverider buoys (X1 to X6) are indicated with triangles. 
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Figure 3. Diagram of a Datawell Directional Waverider buoy and mooring (Datawell bv, 
1995). | 
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SHOWEX: Hurricane Tracks 


Latitude 





Longitude 


Figure 4. Hurricane tracks for Floyd, Gert, Irene and José in 1999. 


39 








































." * FRF Obs a 
— Al): Pasessddncet anaes : alk a H spose = NOAA Buoy 44014 E 30 
OU | | Oo 
© 108 Ceorererarr aren Far as Bee cere eet as @ 0 
é yen On é 
9 se bec cee eee uae agli os... : 9 10 | 
= 0 3 nner = 0 3 
09/15 09/16 09/17 09/18 10/17 10/17 10/18 10/18 10/19 
_~ 360 _~ 360 
7 7 
= 210 = 270 
Cc | a 
= 180 2 180 
: : 
oO Oy ra ai 
a Oo a 
= 0 = 90 — 
= 09/15 09/16 09/17 09/18 = 10/17 10/17 10/18 10/18 10/19 
10 10 
Fl scart eed eee easter eos sateen: 7 cag Pe earercnn eo Rear ore eta) mn (ary Saree mr Breer 
Be Bowe AA hl pet ate © He psec eran 
r rc 
95 ioe ok re On eS FY ECE” 0 SURRTEESS (OTT PERREEESTEST SORTED 
0 : 5 0 : : 
09/15 09/16 09/17 09/18 10/17 10/17 10/18 10/18 10/19 
Date Date 


Figure 5. Time series of wind and wave observations during a) Hurricane Floyd: 
September 15-18, 1999, and b) Hurricane Irene: October 17-19, 1999. Wind speed 
observations from the FRF station and NOAA buoy 44014 (top), wind direction (middle) 
and significant wave height, H, at the offshore buoy, X6 (bottom). The vertical lines 
indicate the time period analyzed for each event. 
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Figure 6. Time series of wind and wave observations during a) a storm on November 2- 
4, 1999, and b) a nor’easter storm on November 10-13, 1999. Wind speed observations 
from the FRF station and NOAA Buoy 44014 (top), wind direction (middle), and 
significant wave height, H;, at the offshore buoy, X6 (bottom). The vee lines indicate 
the time period analyzed for each event. 
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Figure 7. a) The observed ratio 
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Figure 8. Scatter plots of the second-order (a2, b2) vs first-order (al, b1) approximations 
of a) bulk mean propagation direction, Omean, and b) bulk directional spread, op. 
Estimates are compared for the 12 case studies listed in Table 2. | 
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Figure 10. a) Bulk mean direction, 8mean, observed in easterly swell cases vs longitude. 
b) Bulk directional spread, o9, for the same cases vs longitude. 
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Figure 11. a) Hg, b) @mean, and c) Gg vs longitude for large oblique angle swell cases. 
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Figure 13. a) Hs, b) Omean, and c) 6g vs longitude for nor’easter swell case. 
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Figure 14. Frequency and direction spectra of all 6 buoys (X1-X6) for Hurricane Floyd. 
a) The energy density spectra, E(f); b) mean propagation direction, Omean (/); and 
c) the directional spread, og (/). 
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Figure 15. Frequency and direction spectra of all 6 buoys (X1-X6) for Hurricane Irene. 
a) The energy density spectra, E(f); b) mean propagation direction, 8mean (); and c) the 
directional spread, oe (/). 
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Figure 16. Frequency and direction spectra of all 6 buoys (X1-X6) for the November 2 
storm. a) The energy density spectra, E(f); b) mean propagation direction, Omean (/); 
and c) the directional spread, 9 (/). | 
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Figure 17. Frequency and direction spectra of all 6 buoys (X1-X6) for the nor’ easter 

case. a) The energy density spectra, E(f); b) mean propagation direction, Omean (/); and 

c) the directional spread, oo (/). 
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